      SUBROUTINE SETUP
      RETURN
      END
C   HS33
      SUBROUTINE EF(X,FX)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),FX
      ICF=ICF+1
      FX=(X(1)-1.D0)*(X(1)-2.D0)*(X(1)-3.D0)+X(3)
      RETURN
      END
      SUBROUTINE EGRADF(X,GRADF)
      INCLUDE 'FUCO.INC'
      DOUBLE PRECISION X(NX),GRADF(NX)
      ICGF=ICGF+1
      GRADF(1)= 3.0D0* X(1)**2-12.D0*X(1)+11.D0
      GRADF(2)= 0.D0
      GRADF(3)= 1.D0
      RETURN
      END
      SUBROUTINE EH(I,X,HXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),HXI
      CRES(I)=CRES(I)+1
      RETURN
      END
      SUBROUTINE EGRADH(I,X,GRADHI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GRADHI(NX)
      CGRES(I)=CGRES(I)+1
      RETURN
      END
      SUBROUTINE EG(I,X,GXI)
      INCLUDE 'FUCO.INC'
      INTEGER I
      DOUBLE PRECISION X(NX),GXI
      IF ( GUNIT(1,I+NH) .EQ. -1 )CRES(I+NH)=CRES(I+NH)+1      
      IF(I .GT. 3)       GOTO 100
      GXI=X(I)
      RETURN
  100 CONTINUE
      IF(I-5)    200,300,400
 200  CONTINUE
      GXI=5.D0-X(3)
      RETURN
 300  CONTINUE
      GXI=X(3)**2-X(2)**2-X(1)**2
      RETURN
 400  CONTINUE
      GXI=X(1)**2+X(2)**2+X(3)**2-4.D0
      RETURN
      END
      SUBROUTINE EGRADG(I,X,GRADGI)
      INCLUDE 'FUCO.INC'
      INTEGER I,J
      DOUBLE PRECISION X(NX) ,GRADGI(NX)
      IF ( GUNIT(1,I+NH) .EQ. 1 )RETURN
      CGRES(I+NH)=CGRES(I+NH)+1
      IF(I .GT. 4)       GOTO 200
      DO      10      J=1,3
      GRADGI(J)=0.D0
  10  CONTINUE
      IF(I .EQ. 4)      GOTO 100
      GRADGI(I)=1.D0
      RETURN
  100 CONTINUE
      GRADGI(3)=-1.D0
      RETURN
  200 CONTINUE
      IF(I .GT. 5)      GOTO 300
      GRADGI(1)=-2.D0*X(1)
      GRADGI(2)=-2.D0*X(2)
      GRADGI(3)= 2.D0*X(3)
      RETURN
  300 CONTINUE
      DO      400      J=1,3
      GRADGI(J)=2.D0*X(J)
  400 CONTINUE
      RETURN
      END
      BLOCK DATA
      INCLUDE 'BLOCO.INC'
      INTEGER I,J
      DATA NAME/'HS33'/
      DATA X(1),X(2),X(3)/2* 0.D0,3.D0/
      DATA N/3/ , NH/0/ , NG/ 6/
      DATA DEL0/1.00D0/ ,TAU0/1.D0/ ,TAU/.1D0/
      DATA ((GUNIT(I,J),I=1,3),J=0,6)
     F    /-1,0,0,1,1,1,1,2,1,1,3,1,1,3,-1,-1,0,0,-1,0,0/
      END
